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STELLAR ACTIVITY AND TRANSITS 


G. Bruno! and M. Deleuil? 


Abstract. From an observational standpoint, stellar activity poses a 
critical challenge to exoplanet science, as it inhibits the detection of 
planets and the precise measurement of their parameters. Radial ve- 
locity and transit searches revealed a significant fraction of exoplanet 
hosts is active, and showed the need to fully understand the differ- 
ent facets of stellar activity and its impact on observables. Moreover, 
the activity correction is of prime importance for the detection and 
characterisation of Earth analogues. We present a review of the ef- 
fects that stellar activity features such as starspots, faculae, and stellar 
granulation have on photometric and low-resolution spectroscopic ob- 
servations of exoplanets, and discuss the main aspects of the techniques 
which were developed to reduce their impact. 


1 Introduction 


Our understanding of planetary systems is rooted in the knowledge of host stars’ 
parameters and characteristics. Starting from the protoplanetary disc phase, stel- 
lar properties affect the types of planets that will form, as well as their composition 
(e.g. Madhusudhan et al. |2017; Mordasini [2018). As systems age, dynamical and 
magnetic interactions, together with stellar irradiation, are thought to importantly 
affect the evolution of planetary bodies (e.g. Lanza |2018). 

The precision obtained on stellar parameters also has a direct impact on the 
measure of planetary parameters, namely their masses and radii, which are derived 
relatively to the same parameters of their hosts (e.g. Santos and Buchhave |2018). 
In turn, the masses and radii of exoplanets are necessary to calculate their mean 
densities, model their interiors, interpret their transmission spectra and compute 
atmospheric pressure scale heights (e.g. Seager et al. Guillot and Gautier 


Dorn et al. Brugger et al. [2017). 
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While from one side stellar parameters must be determined with high accuracy 
and precision, stochastic stellar brightness variations also hamper the precision 
on measured quantities and prevent a correct determination of planetary param- 
eters and properties. This might result in biases in planetary radii, masses and 
atmospheric composition. 

In this review, the issue of stellar activity from the point of view of photometric 
and low-resolution transit observations is presented and discussed. In Section[2| we 
introduce some of the main stellar activity indicators which are used for exoplanet 
surveys. In Section [3] we show some benchmark cases of active stars observed in 
the context of exoplanet transit surveys, as well as the main techniques used to 
constrain the properties of starspots and faculae from observations. Sections [4] 
and [6]are dedicated to the effect of stellar activity features on transits observed in 
photometry and low-resolution spectroscopy, and the impact of stellar granulation 
is the focus of Section [7| Section [8| presents some conclusions and perspectives 
which can be drawn from what we have learnt so far. 


2 Solar and stellar activity observations 


The first observations of stellar brightness variations come from the best-studied 
star we know, the Sun. Hints of the presence of darker regions on the solar disc, 
which we now call “sunspots”, were recorded in China around 800 BC. In the west- 
ern world, we can find the oldest records in the drawings in John of Worcester’s 
Chronicles, dating 1128. Thanks to the pioneering use of the telescope for astro- 
nomical observations, Galileo Galilei spotted inhomogeneities on the solar disc and 
recorded them on a plate dating 1612, as shown in Figure [I] The observation of 
stellar activity phenomena in the visible part of the spectrum, which were later 
extended to a broader range of wavelengths, enabled also the first detection of a 
solar flare by R. C. Carrington in 1859. 

Pioneering spectroscopic observations of Call H&K line emission at 393 and 397 
nm in solar-type dwarfs were obtained by Wilson and Noyes et al. 
Time variations were shown for other spectroscopic lines (e.g. Livingston et al.|2007), 
such as Ha. Figure |2| presents Call H&K lines of a non-active and an active 
solar-type star, whose spectral line core is affected by strong emission, seen as an 
indicator of stellar magnetic activity. These activity-related features revealed a 
dependence of the chromospheric activity level on the spectral type and the rota- 
tion period of the star (Wright et al. (2004). Main-sequence cool dwarfs have their 
rotation braked due to the loss of angular momentum through a magnetised stellar 
wind. Hence, rotational velocity generally decreases with age, as well as activity 
(Gallet and Bouvier Salabert et al. |2016), even if cases were detected of 
possible weakened magnetic braking and anomalously rapid rotation at old ages 
(van Saders et al. (2016). All in all, chromospheric activity is often used as an 
indicator of age (e.g. Pace and Pasquini Metcalfe and van Saders{2017). 

Photometry also provides indications of time-dependent, stellar brightness vari- 
ations due to the passage of dark starspots and bright faculae across stellar discs. 
The amplitude of such variations clearly correlates with spectroscopic indicators 
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Fig. 1. Illustration from Galileo’s [storia e dimostrazioni intorno alle macchie solari e 
loro accidenti (“History and Demonstrations Concerning Sunspots and Their Properties” , 
or “Letters on Sunspots”), 1613. Downloaded from Encyclopædia Britannica (Pallardy 
https://www.britannica.com/science/sunspot). 
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Fig. 2. The Call H and K lines in the SOPHIE spectrum of a non-active (top) and 


active (bottom) star. From Boisse et al. (2010), A&A, 523, A88, p. 9 (their Figure C.1), 
reproduced with permission ©ESO. 


(e.g. Baliunas and Soon}1995). Additionally, observations provided evidence that 
the level of activity decreases with stellar age, corresponding to a transition from 


spot-dominated to facula-dominated stars (e.g. Lockwood et al. |2007; Preminger, 
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Telescope Precision Time sampling mag limits Sample size | Reference 
[stars] 
SOHO/VIRGO < 10 ppm/min 60 s ~ The Sun Frohlich et al. (1995), 
photon noise Frohlich et al. (1997), 
Jiménez, Roca Cortés, and Jiménez-Reyes (20 
MOST ~ 100 ppm/hr 1— 60s 04<V <6 5000 
CoRoT 700 ppm/hr 32 or 512 s 11.5< V <16 170000 
Kepler/K2 80 ppm/hr 59 s or 29 min T< Kp<17 ~ 170000 
TESS 60 ppm/hr 120 s or 30 min 4<V <12 > 200000 
CHEOPS 50 ppm/hr 1.05 — 60 s 6<V <13 > 600 
PLATO 27 ppm/hr (goal) 2.5 or 25 s 8<V<16 > 108 
Ariel 10-100 ppm/hr (goal) z bright stars ~ 1000 


Table 1. Space missions dedicated to ultra-high precision stellar photometry and spec- 
trophotometry. 


Chapman, and Cookson {2011} Shapiro et al. Radick et al. Reinhold 
et al. (2019). In this regard, our Sun belongs to the category of weakly-active stars 
(Lockwood, Skiff, and Radick |1997). Photometric modulations are observed on 
various time scales, from days or tens of days — indicating starspots and faculae 
brought about by stellar rotation — to months or years — revealing the evolution 
of the brightness dishomogeneities pattern on the stellar disc, activity cycles, and 
variations in the stellar disc coverage of activity features, or “filling factor” (e.g. 
Nielsen et al. (2019). This information enables placing constraints on stellar ro- 
tation rates, and on the distribution and evolution of activity features, and is 
therefore particularly helpful to understand the underlying dynamo mechanism 


across different spectral types (e.g. Hall}1991} Lanza et al.|1998} Berdyugina|2004 
2005} Strassmeier |2005) |2009; Bradshaw and Hartigan |2014} See et al./2016} Brun 
and Browning |2017; Herbst et al. |2021). 


In order to monitor daily photometric variations, ground-based observations 
are however not sufficient. Day-night alternation forces continuous interruptions, 
with adverse weather and difficulty to obtain continuous, long-term telescope time 
worsening the issue. In order to assess the impact of stellar brightness variations 
on transit observations, long-term uninterrupted monitoring is necessary. Hence, 
both exoplanet and stellar activity science greatly benefited from the launch of 
space telescopes dedicated to ultra-high precision and long-duration photometric 
observations, and better characterised the time scales concerned by stellar activity 


features (e.g. Ballot et al. |2011; Bradshaw and Hartigan |2014). Table [1] gives a 
summary of these space-borne telescopes. 


3 High precision photometry 


Space missions provided ultra-precise long-duration observations of hundreds of 
thousands of star, covering nearly all spectral types. However, even for the more 
than four year-long duration of the Kepler observations, the duration of these 
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Fig. 3. Examples of solar light curves (upper left) compared with sample Kepler light 
curves for solar-like (upper right), quiet (lower left), and spotted (lower right) stars. 
The normalised light curves are offset with respect to each other for clarity. Solar light 
curves are coloured in red, and Kepler light curves in black. In each panel, a solar light 
curve is shown for comparison. An offset was applied to each light curve. From Basri 
et al. (2010), ApJ, 769, 37, pp. L155-L159 (their Figure 2). Reproduced by permission 
of the AAS. 


missions is too short to cover the entirety of a stellar cycle (which, for the Sun, 
lasts about 11 years). Despite the limited temporal coverage, it was possible to 
observe a wide variety of brightness variation patterns and to connect them to 
stellar types. Figure |3] juxtaposes some light curves of the Sun to those of some 
Kepler stars grouped by activity level. Comparative studies of this kind placed the 
Sun among weakly-active stars, or in a transition phase between active and non- 
active stars (Basri et al. Basri, Walkowicz, and Reiners Bohm-Vitense 
Brandenburg, Mathur, and Metcalfe Strugarek et al./2017). Moreover, 
observations on as much as ~ 10° stars confirmed expectations that cool, mainly 
convective stars exhibit larger photometric variations (e.g. García et al. [2014). 
This flux modulations observed in long-term and continuous light curves can 


reveal surface feature topologies and properties. As an example, Figure [4] shows 
the light curve of CoRoT-2, hosting one of the first exoplanets discovered by the 
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Fig. 4. The light curve of CoRoT-2, where the transits of the planet are evident as 
vertical flux drops. In cases like this, the out-of-transit stellar flux cannot be determined 
with high confidence. Credit: Alonso et al. (2008), A&A 482, L21-L24, p. 3 (their Figure 
1), reproduced with permission QESO. 


CoRoT mission. From photometric amplitude variations, the size and brightness 
contrast (i.e., temperature) of starspots and faculae can be extracted. Deviations 
from a simple sinusoidal shape, which would be an indication of a single activity 
feature, increase as the number of stellar activity features increases and as their 
shape gets more complex. We have to remember, however, that we are not able 
to resolve the stellar disc, so that we cannot distinguish between large stellar 
activity features and agglomerations of smaller features. In addition, the complex 
inter-dependence of rotation, differential rotation, and the magnetic field-induced 
features that evolve in time, produces degenerate cases that hamper the inversion 
(e.g. Mosser et al. (2009). Despite this, long-term data series that allow a time 
coverage of a few stellar rotations have triggered new developments, bringing new 
insights into the configuration and evolution of starspot and faculae, but also key 
constraints on the stellar rotation period and differential rotation. 


3.1 Stellar rotation periods 


Spots or inhomogeneities on the stellar surface are dragged away by the rotation 
of the star, resulting in a modulation of its brightness over days or tens of days. 
This modulation can be used to infer its rotation period, and was taken advantage 
of by photometric surveys which made large samples analyses possible (e.g. Basri 
et al. {2010} Affer et al. (2012). 

By applying the autocorrelation function (proven to be more robust than the 
Lomb-Scargle periodogram for this kind of problem), McQuillan, Mazeh, and 
Aigrain were able to measure the rotation period of more than 34 000 stars, 
over a sample of more than 133 000 stars. Among other conclusions, they reported 
an increase of the rotation period with decreasing stellar mass and temperature, 
broadly consistent with the stellar spin-down law. Their findings, later expanded 
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Fig. 5. Left: Geometric representation of a planet (black) transiting in front of a group 
of active regions (red). The axes denote the position on the star, normalised to the stellar 
radius, and are centred on the centre of the stellar disc. Right: Corresponding normalised 
light curve. From Montalto et al. (2014), MNRAS 444, 1721-1728 (their Figure 2). 


to an even larger sample (e.g. Lanzafame et al. Savanov and Dmitrienko 
Santos et al./2019) also suggest that cooler stars are more evenly covered in 
brightness inhomogeneities such as starspots — a factor that is particularly crucial 
to discern the effect of stellar activity on exoplanet transits. This aspect will be 
better explored in Section [5] 


3.2 Stellar surface reconstruction 


Knowledge on the distribution of stellar surface inhomogeneities can provide very 
valuable insights on the amount by which transit parameters are affected. Luckily, 
photometric light curves encode information on the distribution and temperature 
of such features. The complex problem of light curve inversion has been tackled in 
a number of ways, since much earlier than it became relevant to exoplanet science 
(e.g. Budding Eaton and Hall [1979} Vogt Rodonò et al. |1986). Here 
we report some of the main approaches that have been used in the literature. 


e Analytic models describe the stellar disc as the visible part of a sphere emit- 
ting a total unit flux, from which the contribution of darker (or brighter) 
regions with circular shape is subtracted (or o to a the light 
curve wo E Dorren [1987} AA 4; Kipping |2 Montalto et 
al. 4; Béky, Kipping, and Holman |2014} Ikuta et al. oer Usually, 
m inhomogeneities are nan por according to thio projection at a 
given stellar longitude and latitude. Other parameters can be tuned, such 
as the inclination of the stellar rotation axis, the stellar rotation period, 
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or the stellar limb darkening coefficients. Such a “forward” model is then 
implemented in a Markov Chain Monte Carlo (MCMC) algorithm (or its 
variations, such as Nested Sampling) to find the posterior probability dis- 
tributions for the model parameters (e.g Kipping Bruno et al. 
Ikuta et al. 

An analytic model requires only a fraction of a second to be calculated 
and, in principle, there is no limit to the number of stellar activity features 
that can be modelled. Clearly, the more complex the model, the longer the 
computation time required to reach convergence. It is then often desirable 
to adopt some simplifications, as degenerate cases among parameters might 
hamper the convergence of the chains. In particular, not all models allow the 
modelling of planet-starspot occultations, which provide stong constraints 
on the active region distribution along the transit chord (see Section (4.2). 
Dedicated modelling approaches were developed for this specific problem 
(Tregloan-Reed, Southworth, and Tappert |2013; Béky, Kipping, and Hol- 
man |2014; Montalto et al. Scandariato et al.|2017} Juvan et al. (2018). 
As an example, Figure |5| represents the modelling of a planet occulting a 
particularly complex combination of bordering active regions. 

However, observations provide no constraint on some of the parameters. To 
overcome this issue, the temperature of the features can be fixed by adopt- 
ing results from the literature (e.g. Berdyugina [2005), or the stellar rotation 
axis can be fixed when there is no constraint on such parameter. Moreover, 
a large number of starspots or faculae might overcomplicate the model, and 
prevent a proper sampling of the parameter space by the MCMC. 


e Alternatively, the stellar visible disc can be modelled with a grid of hundreds 

of elements, and the contribution of each cell to the total flux can then 
be evaluated. This allows the modelling of complex shapes for starspots 
and faculae, but is computationally very expensive (e.g. Boisse, Bonfils, and 
Santos {2012 Oshagh et al. [2013) P] 
The activity feature distribution can be retrieved in different ways. One 
is Maximum Entropy Regularisation, where the brightness contrast (or the 
effective temperature) of starspots and faculae for each cell of the grid is 
fixed, and a functional that models the light curve given the filling factor f 
of the cells (i.e., their fraction covered by activity features) is minimised: 


© = X? (f) — AS(f), (3.1) 


1. Among publicly available codes: macula (https://github.com/davidkipping/MACULA 
Kipping |2012), KSint (http://eduscisoft.com/KSINT/| Montalto et al. |2014), Prism (https: 
//github.com/JTregloanReed/PRISM_GEMC, Tregloan-Reed, Southworth, and Tappert 
and TOSC Scandariato et al. |2017). 

2. In particular, the SOAP 2.0 code (Boisse, Bonfils, and Santos Oshagh et al. 2013) 
allows such a modelling, and is discussed in other proceedings of the EES2019. The code was 
presented by M. Oshagh at the school (https://ees2019.sciencesconf.org/resource/page/id/9) and 


can be downloaded at http://www.astro.up.pt/resources/soap-t/ (installation instructions at 
https: //ees2019.sciencesconf.org/resource/page/id/8). 
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Fig. 6. Reconstruction of the active region longitudinal distribution (x-axis) using hun- 
dreds of days of Kepler-17 photometric data (time on the y-axis). Modified from Lanza 
et al. (2019), A&A 626, A38 (their Figure 9). The colours, from blue to orange, indicate 
an increasing filling factor, and identify active longitudes. Overplotted in white, the 
active region distribution identified by Valio et al. using occulted starspots (see 
Section [4.2). Reproduced with permission ©)ESO. 


where y? is the usual chi-square, S is the entropy functional (a function 
describing the complexity of the model, so that more complex models are 
penalised) and A is a Lagrange multiplier. Despite the large number of 
free parameters, stable solutions can be achieved with this method, which 
retrieved active longitudes from some CoRoT and Kepler light curves of 
exoplanet-host stars (e.g. Lanza et al. |2009; Lanza et al. (2019). Figure [6] 
shows the reconstruction of the active region distribution using the light 
curve of Kepler-17 (Lanza et al. [2019). Another option is to use a relatively 
small number of cells to describe the stellar disc: this allows one to adopt 
variations of the x? statistics, or even MCMC algorithms to achieve the best 


fit (e.g. Wolter et al.|2009} Huber et al. |2010; Tregloan-Reed et al. |2015). 


3.3 Stellar differential rotation and starspot evolution 


In order to reconstruct the properties of the stellar dynamo and the generation 
of magnetic fields in the stellar convection zone, efforts have been directed to 
collect information on the variation of the stellar rotation period at different stellar 
latitudes (e.g. Knobloch and Spruit Brandenburg and Subramanian 
Balona and Abedigamba |2016). This variation can be expressed as 


Q = Q — AQ sin? Y, (3.2) 


where ~ is the stellar latitude, Q the rotation rate at the stellar equator and AQ 
the variation of the rotation rate between the stellar pole and the equator. The 
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Fig. 7. Reconstructed light curve map of GJ 1243, from Davenport, Hebb, and Hawley 
(2015), ApJ 806, no. 2, 212 (their Figure 5). The x-axis corresponds to the observation 
date, and the y-axis to the stellar rotation phase. The two retrieved starspots are marked 
in orange and purple. The orange pattern represents the higher-latitude, more stable 
starspot, and the purple one indicates the lower latitude feature migrating in longitude 
because of differential rotation. Reproduced by permission of the AAS. 


“shear” a = Qo/AQ is one way to express the deviation of stellar rotation from 
the one of a rigid body, and is equal to 0.2 for the Sun (AQọ = 0.055 rad day” *). 

A variety of methods exist to estimate differential rotation from photometry 
and spectroscopy (e.g. Berdyugina Lanza, Das Chagas, and De Medeiros 
Nielsen et al. Santos et al.|2017). The advent of space-borne high 
precision, long-term photometry added a powerful resource, both for the quantity 
of continuous data (e.g. Lanza|2014} Garcia et al.|2014), allowing to precisely track 
stellar rotation periods as a function of stellar latitude, and to asteroseismology 
(Benomar et al. Bazot et al. |2019). Here we discuss the example of Dav- 
enport, Hebb, and Hawley (2015), who applied a starspot modelling technique to 
reconstruct the longitudes and radii of starspots on the Kepler light curve of the 
M dwarf GJ 1243. Using 5 days-long sliding windows and an MCMC sampler, 
they derived longitudes and radii of two starspot features: a high-latitude feature 
for over 6 years of observation, and a low-latitude feature evolving in both size 
and longitude over hundreds of days. For this second feature, they measured a 
differential rotation rate AQ = 0.012 + 0.002 rad day’ and a shear a = 0.00114, 
one of the lowest values measured for cool stars. The phased light curve map is 
illustrated in Figure [7] 

However, light curve modulations are degenerate with respect to the latitude 
of stellar activity features, so that these parameters were fixed during the mod- 
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elling. As discussed in Section [4.2] the occultation of starspots and faculae during 
transits allows setting constraints on their latitudes, and hence provides significant 
information for better understanding the dynamo of the host star. Starspot and, 
for the first time, facula occultations were used to derive the differential rotation 
rate of a young solar-type star, Kepler-71, by Zaleski et al. (2019). The result 
of AQ < 0.005 rad day~!, much lower than that of the Sun’s, describes a rigidly 
rotating star which is set apart from solar-type stars at a similar evolutionary 
state. 

Another major limitation when trying to estimate the configuration of activity 
features is the description of their evolution, as most of the usual models are 
stationary. One way to address this issue is to split the light curve in chunks and 
to model them with different distributions of starspots and faculae (Silva-Valio et 
al. Lanza et al.|2011} Davenport, Hebb, and Hawley |2015). Indeed, stellar 
activity features do evolve in time, and characteristic laws can be used to describe 
their growth from first appearance to maximum size, permanence at maximum 
size and decay until disappearance (e.g. Kipping |2012} Bradshaw and Hartigan 
van Driel-Gesztelyi and Green |2015} Savanov Namekata et al. |2019). 
With a starspot and facula model including feature evolution, longer segments of 
the light curve can be modelled and the lifetime of individual features (or groups) 


can be inferred (e.g. Mosser et al. |2009; Bruno et al. |2016). 


4 Planetary transits and stellar activity features in photometry 


Stellar activity directly affects planetary transits in different ways, as spots and 
faculae leave their imprint on both their depth and shape. This impacts the 
observables used to derive the radius of an exoplanet, and as such it is of crucial 
importance to disentangle them from any kind of contamination. 


4.1 Effect of non-occulted activity features on the transit depth 


Variations in the apparent transit depth are caused by changes in the out-of-transit 
stellar brightness. If Fy, is the out-of-transit stellar flux and Fin the in-transit 
stellar flux (both at the centre of the transit), the transit depth is defined as 


Fou, — F ? 
pa ont ni -(}) l (4.1) 
Ry 


where rp and R, are the planetary and the stellar radius, respectively. If the star 
is covered by activity features, Fout — Fin remains the same, but Fout changes with 
time: in particular, if it decreases (increases) because of starspots (faculae), the 
transit depth will increase (decrease). Brightness variations are partly accounted 
for by fitting the out-of-transit flux in proximity of the transit edges with a low- 
order polynomial, and by computing its value along the transit. The relative effect 
is more remarkable the larger the brightness modulations are, and the smaller the 
planet is. Variations can span from a few 10? ppm to a few percent (Czesla 
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et al. 2009), with considerable consequences on the inferred planetary properties 
(such as its internal structure and composition, e.g. Wagner et al. (2011). Czesla 
et al. (2009) suggested that a different quantity than the transit depth should be 
adopted in order to minimise the effect of stellar brightness variations. In their 
notation, f = Fout — Fin, n = Fout, and the proposed metric is 


fin 
P 


where p is the value of the stellar photometric flux unaffected by starspots or 
faculae. The problem with this expression, however, is that the contribution of 
faculae to the stellar flux is not known a priori, so that the value of p cannot be 
deduced from observations nor theory. 

Another way to estimate the out-of-transit stellar flux is via the use of starspot 
and faculae modelling. If observations covering a few stellar rotations are avail- 
able, the distribution of starspots and faculae, as well as their brightness contrast 
(i.e., effective temperature) can be reconstructed up to some level of confidence, 
and Fout at any intermediate time can be inferred. The precision on this value 
depends on the uncertainty on the retrieved activity feature distribution and its 
degenerate solutions. One can remove the transits from the light curve, describe 
the out-of-transit flux with one of the models outlined in Section B.2] and use this 
value for the in-transit flux at any given time. Alternatively, models enabling si- 
multaneous simulations of both activity features and transits can be fitted to the 
entire light curve. In this way, constraints on both the transit and the activity fea- 
ture parameters can be obtained, and the need for transit normalisation avoided 
(Bruno et al.|2016). The price for this is a considerable increase in computational 
and model selection effort for the transit fit. 

A third option, which has recently gained popularity, is to apply a Gaussian 
process (GP) regression to the whole light curve altogether, setting the transit as 
the mean function of the process (Aigrain, Parviainen, and Pope {2016} Foreman- 
Mackey et al. Barros et al. With a GP, the stellar signal is modelled 
as correlated noise, without the need to describe it explicitly. However, this re- 
quires the choice of the functional form for the correlation among observations, 
i.e. a “kernel”. Some specific kernels are often adopted to model the variations 
induced by starspots and faculae, such as combinations of the exponential-squared 
and exponential-sine-squared kernels (Haywood et al. (2014). The GP can be cali- 
brated so that the kernel parameters (or “hyperparameters” ) are used to perform 
posterior inference on the size, position, and contrast of active regions, as well 
as their evolution time scales (Luger et al. Luger, Foreman-Mackey, and 
Hedges |2021). Convenient open-source software for the GP implementations here 
described is available onlin] and can be readily included into an MCMC simula- 
tion. GPs are particularly powerful thanks to their flexibility, but they are also 


D= 


+1, (4.2) 


3. A presentation dedicated to GP regression was given by S. Aigrain at the EES2019 


//ees2019.sciencesconf.org/resource/page/id/9 
4. Ambikasaran et al. |2015)), 


or example, ub. com / dfm / george 
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very demanding in terms of computation resources: given N data points, this is 
O(N?) for most kernels, and O(N) for a specific kind of kernels applied to one- 
dimensional data sets (e.g. Ambikasaran et al. 2015} Foreman-Mackey et al. [2017). 

Which approach is the most advantageous depends on the case at hand. The 
simplest model might be the most useful when only a few transits are available 
(even only one, as it is sometimes the case). Analytic modelling of starspots 
and faculae is recommended to obtain a physical representation of the starspot 
configuration and temperature, given an uninterrupted light curve, and it can 
provide with a more robust assessment of the out-of-transit flux. GPs, on the other 
hand, avoid the need of precise parametrisation of the geometry of the problem 
and can be coupled with models of instrumental systematics with a relatively low 
number of free parameters. None of these methods is completely exhaustive, and 
the next Section focuses on how occultations of activity features can help break 
some parameter degeneracy. 


4.2 Effect of occulted activity features 


When starspots or faculae are occulted by a planet during a transit, they produce 
typical “bumps” on the transit profile. They affect transit depth measurements 
in an opposite way to non-occulted ones: if the planet passes in front of a dark 
starspot, the apparent transit depth will temporarily decrease, and the opposite 
will happen in case a bright facula is occulted (Figure[gp. From the Sun, we know 
that faculae are often associated to sunspots (Kiepenheuer|1953). They are mostly 
visible when they are close to the stellar limb, and are therefore described by a 
wavelength-dependent limb brightening law (e.g. Spruit {1976} Shapiro, Peter, and 
Solanki {2019}. 

If the temperature contrast of the activity features with the average stellar 
photosphere is large enough, given the noise level of the data, their occultations 
can be easily identified in the transit profile. In principle, they can therefore be 
removed from the transit profile, which can then be fitted with a standard transit 
model. However, if their signal cannot be clearly singled out, or during a given 
single transit the planet crosses multiple stellar activity features, removing the 
imprint of the occultation might be problematic. The case of undetected occulted 
features was discussed by e.g. Ballerini et al. (2012), while the problems of multiple 
starspot crossings in the same transit were explored by e.g. Czesla et al. 
and Silva-Valio et al. (2010) on the benchmark case of the hot Jupiter CoRoT-2b. 
These latter authors found that, under the assumption that all crossed features 
are dark starspots, the transit radius of CoRoT-2b might be overestimated by up 
to 3% if the crossing effect is not taken into account. However, our knowledge 
of starspots and faculae in other stars is still not complete enough to exclude the 
contribution of faculae, which could average out this effect, or even increase the 


celerite (https://github.com/dfm/celerite 


sors in scikit-learn (https://scikit-learn.org 


Foreman-Mackey et al. (2017) and the GP regres- 
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Fig. 8. Effect of the planet crossing of a spot (left) and a facula (right) on the transit 
profile. From Espinoza et al. (2019), MNRAS 482, no. 2 (their Figure 7). 


apparent transit depth in individual transits. Indeed, Bruno et al. (2016) found 
indications that some of the deepest transits of CoRoT-2b, previously considered as 
the closest ones to the “unperturbed” transit profile, might be affected by faculae. 

A very interesting advantage of starspot and facula crossings is that their po- 
sition on the stellar disc (assuming the stellar rotation axis is known) can be 
determined with a much larger precision than in the case of non-occulted features. 
Provided a stellar inclination angle of 90° with respect to the line of sight, the 
feature latitude ¢ can in fact be constrained within the stellar latitudes crossed 
by the transit chord: 

Tp 
o = arctan (=) . (4.3) 
Thanks to this, the latitudinal distribution of the starspots on HAT-P-11 was 
shown to follow a solar-like behaviour (Morris et al. |2017). 

Occultations also provide precise constraints on the longitude of active regions, 
which are found to be consistent with the reconstruction of the stellar surface 
using only unocculted features (Lanza et al. |2019| see Figure (6). In the special 
case of Kepler-7b, Désert et al. were able to observe the recurrent passage 
of the same group of occulted starspots on the transit profile. Together with a 
previous measurement of the stellar obliquity via the Rossiter-McLaughlin effect, 
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and helped by the commensurability of the stellar rotation and planet revolution 
period (P,/Ppianet œ 8), it was possible to follow the starspots moving across 
different transits. Thanks to the long duration of the Kepler-7 light curve, starspot 
groups were observed without interruption for about 100 days, and important 
constraints on the lifetime of these features were obtained. 

A last aspect to underline is the possibility of constraining the size and temper- 
ature of spots and faculae thanks to their occultation. A minimum value on their 
size can be inferred from the duration of the bump they produce on the transit 
profile, and this can be scaled to a feature size if the transit duration and stellar 
radius are known. With a simple geometric model that takes the activity feature 
brightness contrast with the stellar surface as a free parameter, one can derive the 
temperature contrast of the active feature to the stellar photosphere by inverting 
the ratio between black body laws (Silva [2003): 


_ exp (hv/KpTer) — 1 
~ exp (hv/KpTteat) — 1’ 


(4.4) 


where p is the brightness contrast with the stellar photosphere, h is Planck’s 
constant, Kp is Boltzmann’s constant, v is the frequency of observation and Tog 
and Treat are the stellar and the activity feature effective temperature, respectively. 
This method was adopted in several cases (e.g. Silva-Valio et al. Mancini 
et al. (2017). 

Clearly, the black body law is only a useful approximation as stars, starspots, 
and faculae are not black bodies at different temperatures, and each one has its 
own spectrum and specific absorption lines. Section |5| contains a discussion on 
how transmission spectroscopy can help us probe not only planetary spectra, but 
also those of occulted stellar activity features. 


5 Starspots and faculae in transmission spectroscopy 


Transiting exoplanets are key targets to probe the physical structures and chemi- 
cal compositions of their atmospheres. As their upper atmosphere absorbs and 
scatters starlight, transmission spectroscopy can be used to infer the physical 
conditions within their atmospheres by observing the transit depth at different 
wavelengths (Charbonneau et al. (2002). Here, also, the impact of starspots and 
faculae cannot be neglected because, as discussed in Section they increase 
and decrease the measured transit depth, respectively. The variations are not 
only time-dependent, but also wavelength-dependent: according to Wien’s law, 
the peak emission from a black body moves to longer wavelengths as the black 
body gets cooler, and there is more sensitivity close to the emission peak. Hence, 
the effect of starspots and faculae is stronger in the visible, and less prominent 
(but still important) in the infrared (IR). 

As shown in the example of Figure [9] it was found that unocculted starspots 
(faculae) result in a rising (descending) slope in the transit depth from the near-IR 
to the ultraviolet, which could be wrongly interpreted as Rayleigh-scattering dust 
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Fig. 9. The transmission spectrum of HD 189733 b in different bands. In the visible part 
of the spectrum, blue and orange colours represent Rayleigh scattering and unocculted 
starspot models, respectively. From McCullough et al. (2014), ApJ 791, 55 (their Figure 
4). OAAS. Reproduced with permission. 


(or its absence) in the planet’s upper atmosphere (e.g. Pont et al. Pont et 
al. McCullough et al. {2014} Rackham et al.[2017} Alam et al. Sotzen et 
al. |2020). The opposite happens for occulted active regions (Oshagh et al. (2014). 

Among molecular spectral signatures which are nearly routinely detected, the 
water absorption band is the main one observed by Hubble Space Telescope’s Wide 
Field Camera 3 instrument, and is the main indicator of oxygen abundance and 
metallicity in an exoplanet atmosphere. However, in starspots cooler than ~ 3000 
K, water vapor could exist (Wallace et al.|1995) and produce absorption features 
that can mimic planetary water absorption at ~ 1.4 and ~ 2.3 um (e.g. Wakeford 
et al. |2019). 

Because of this, an assessment of the activity level of the exoplanet host is 
necessary before any interpretation of a transmission spectrum is attempted. The 
correction of the activity-induced transit depth variations can be done by a pho- 
tometric monitoring of the host star, with stellar models (such as Castelli and 
Kurucz 1D ATLAS and PHONENIX Husser et al. to account for the 
wavelength-dependent flux correction, in order to properly evaluate the stellar con- 
tribution (e.g. Huitson et al. Alam et al. (2018). In this kind of modelling, 
the size of the activity features was shown to be essential (Rackham, Apai, and 
Giampapa (2018). With the use of stellar synthetic spectra to model both stars 
with different spectral types and their starspots, Rackham, Apai, and Giampapa 
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Fig. 10. The amount of contamination of transmission spectra for planets hosted by 
active M stars, in the case of giant (left) and solar-size starspots (right), according to 
Rackham, Apai, and Giampapa (2018), ApJ 853, 122 (their Figure 7). The molecular 
absorption bands at various wavelengths are reported. Reproduced by permission of the 


AAS. 


found that active regions with angular size comparable to the size of large 
sunspot groups (~ 2°, Mandal et al. |2017) can produce up to tens of percent 
variations in atomic and molecular absorption bands, against a few percent in 
the case of features which are four times as large. This is due to the fact that, 
for a given activity-induced stellar photometric variability amplitude, small ac- 
tive regions require much larger covering fractions to produce such amplitude, and 
produce a larger contamination in the transmission spectrum. Figure [10] presents 
their findings on spectral contamination from starspots on M dwarfs. With the 
use of similar simulations, the same effects were shown to be much smaller for 
FGK stars (Rackham, Apai, and Giampapa [2019). 


The shape of the out-of-transit light curve could be the key to distinguish giant 
from solar-like activity features on planet hosts, and so to identify optimal targets 
for transmission spectroscopy of planets orbiting cool stars. In the case of giant 
starspots, wide variations of the light curve will be observed, with a periodicity 
which depends on the stellar rotation period. In the case of small starspots, es- 
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pecially if these are uniformly distributed on the stellar disc, brightness variations 
will have a much smaller amplitude and no clear periodicity. This might be the 
most frequent case for M stars (see Section [3) where, on average, a small starspot 
going out of view could be compensated by another one coming into view. The 
two scenarios could be discerned thanks to a continuous monitoring of the host 
star in the time windows during which transits are observed. Inspection of the 
light curve for at least two or three rotation periods in different spectral bands 
should enable to distinguish the type of light curve variations, and so to gauge 
the amount of contamination on the transmission spectrum. Remarkably, a time- 
evolving, longitudinal map and effective temperature reconstruction of the active 
regions on WASP-52 was obtained from ~ 600 days of BV RI multiband photom- 
etry, implying a reduction of their effect on transmission spectroscopy by about 
one order of magnitude (Rosich et al. (2020). A complementary approach would 
be to model the temperature and size of active regions using out-of-transit low- 
resolution spectroscopic observations. This technique showed promising results for 
the correction of transmission spectra affected by large-contrast spots, in modelled 
bright (V = 9) stars (Cracchiolo, Micela, and Peres[2021). 

The time dependence of the “filling factor” (i.e., the fraction of the stellar sur- 
face covered by activity features) also operates a variation in the baseline between 
two or more transit observations. The photometric monitoring of the host star’s 
activity level was, in fact, shown to be key to correctly stitch together multi-epoch 
transmission spectroscopy observations, aiming either at extending the wavelength 
coverage or at detecting temporal variations of the planet’s atmosphere. Retrieval 
exercises including the contribution of activity features to the spectrum were also 
shown to provide constraints on both the atmospheric and the starspot param- 
eters (e.g. Barstow et al. Rackham et al. Zellem et al. |2017; Bruno 
et al. (2020). 

Simultaneous transit observations in multiple wavelengths can provide very 
useful information on the physical properties of activity features when they are 
occulted by the planet. In that case, indeed, the relative flux increase (the height 
of the starspot occultation bump) Af measured in different wavelengths \ probes 
specific absorption lines. This makes it possible to place constraints on the effec- 
tive temperature of the activity feature, whose emission is modelled using stellar 
synthetic spectra. Following Sing et al. (2011), 


Af(A) _ 1— Bis) a 
Afo) 1 F= FR 


where FY is the normalised flux measured at wavelength x and effective temper- 
ature y, and ào is a reference wavelength. Figure shows the so-determined 
spectrum of a starspot occulted during a transit of HD 189733b. For the specific 
case of this active K0-dwarf star, Sing et al. derived a spot temperature 
Tspot ~ 4250 + 250 K and spot-to-star temperature differences compatible with 
what is measured for sunspots. In the case of WASP-19, the fit of the occultation 
distortion at varying wavelengths allowed the determination of better than 100 K 
constraints on a bright and a dark active region (Espinoza et al. [2019}. 


(5.1) 
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Fig. 11. Relative height of the spot occultation bump with respect to the same quantity 
at 600 nm, as a function of wavelength (data points), and Kurucz 1D ATLAS models used 
to fit the measured flux rise during a transit of HD 189733b. Measurements (symbols) 
were obtained from Hubble Space Telescope’s STIS G430L and ACS instruments. Starspot 
models range from 4750 to 3500 K in 250 K intervals (blue to orange, respectively), while 
a Teg = 5000 K ATLAS stellar model (Castelli and Kurucz {2003) was used for the host 
star. From Sing et al. (2011), MNRAS 416, 1443-1455 (their Figure 10). 


6 Effects on other transit parameters 


Among the transit parameters which are affected by stellar activity features, limb 
darkening (LD) is one of the most crucial. In turn, centre-to-limb stellar bright- 
ness variations modify both the transit profile and its apparent depth, affecting 
the determination of the exact radius of the planet (Mandel and Agol [2002} Csiz- 
madia et al. Heller [2019). In Figure [12 left panel, the transit profile shape 
variation is shown when the effect of LD is included. In the figure, the transit 
depth increase is dubbed “overshoot”. On the right, the overshoot is analytically 
derived and shown for different spectral types. If stellar activity features affect 
the LD coefficients, the transit depth measure (and therefore the measured radius 
of the planet) for a given stellar type is prone to error. 


Usually, LD coefficients are obtained by fitting stellar centre-to-limb specific 
intensity model profiles with a variety of increasingly complex analytic relation- 
ships. For example, two popular laws are the quadratic (Kopal and the 
non-linear (Claret LD law. Using 0 as the angle between the stellar surface 
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Fig. 12. Left: the effect of stellar limb darkening on the transit shape and depth. The 
increase in transit depth over the planet-to-star surface ratio is here called “overshoot” 
(o). Right: The variation in overshoot (oint, z-axis) for a quadratic limb darkening 
law (coeffcients a and b) and orbit inclination i = 7/2 rad for different stellar types. 
Variations in the limb darkening coefficients due to stellar activity cause variations in 
the overshoot. From Heller (2019) (Figures 1 and 3), A&A 623, A137, pp. 2 and 5, 
reproduced with permission (Q)ESO. 


normal and the line of sight, they are parametrised as 


D (u) = 2 
TD =1- a(l- u) — a2(1—p)?, (6.1) 
and i 
D (u) _ n/2 
a t aiy =), (6.2) 


respectively, where u = cos@ and Iy(j)/I)(1) expresses the ratio between the 
specific intensity at any wavelength À and position u and the centre of the stellar 
disc (corresponding to u = 1). 

The LD coefficients an can be fitted to the stellar models and then used as fixed 
parameters to model the transit profile (Claret [2000). In particular, the non-linear 
law has been shown to produce the most accurate fit to the models so far (e.g. 
Morello et al. and references therein). 

A fit of the LD coefficients to the transit profile is often more desirable than 
fixing them to theoretical values. As illustrated in Figure[13} left panel, this is due 
to the unavoidable uncertainties in the stellar parameters, as well as to differences 
in the stellar models that are used. This is particularly the case when stellar ac- 
tivity features produce time-dependent slight deviations of the stellar parameters, 
that affect the LD coefficients and therefore the transit depth overshoot. When 
the LD coefficients are included in the transit fit as free parameters, however, the 
non-linear law often hampers convergence of the MCMC sampling of the transit 
parameters, and simpler laws need to be used. Starspots were found to potentially 
produce a few percent variations in the quadratic LD coefficients, as well as up 
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Fig. 13. Left: Variation of k = rp/R, from the overshoot given by Sing and 
the Claret and Bloemen (2011) LD coefficients, for a low (continuous line) and grazing 
(dashed line) impact parameter b. Right: variations in the observed LD coefficients for 
varying starspot filling factors for a star with Teg = 5775 K and Tspot = 3775 K. The 
positions of the spots were randomly chosen on the visible hemisphere. Figures from 
Csizmadia et al. (2013), A&A 549, A9, pp. 4 and 8 (their Figures 2 and 5), reproduced 
with permission (QESO. 


to ~ 10% variations in r,/R,, depending on the stellar effective temperature and 
transit impact parameter (Csizmadia et al./2013). The maximum effect was found 
around 7.4 = 4000 K, grazing transits, and for large spot filling factors, as shown 
on the right panel in Figure [I3] Flux-transport models based on the Sun showed 
that the effect of faculae can also be non-negligible, and that deriving LD curves 
from the transit fits with minimum residuals induces wavelength-dependent biases 
due to both spots and facuale (Schrijver [2020). 

Because of the limited accuracy of the quadratic LD law, some variations were 
proposed, as well as efficient calculation methods and proper priors on the coeffi- 
cients to include in the transit fit (e.g. Kipping {2013} Espinoza and Jordan {2016} 
Morello et al. Claret Maxted|2018). More accurate three-dimensional 
stellar models also help us better understand the impact of granulation and the 
behaviour of LD with fundamental parameters (e.g. Magic et al. (2015). Due to 
variations on individual stars and the time dependence of the activity pattern on 
the stellar photosphere, however, analysing exoplanet transits still requires a com- 
promise between accuracy and computation efficiency when adopting an LD law 
to implement in an MCMC. 

As previously highlighted, LD is a wavelength-dependent effect: it is stronger in 
the visible, but its effects are also measured in the IR. It has therefore to be taken 
into account both in transit photometry and spectrophotometry. Observations 
with the Hubble Space Telescope deserve a particular mention: objects which do 
not lie in the continuous viewing zone of HST cannot generally be observed for 
a full continuous transit. Sometimes, the transit edges are sacrificed in order to 
obtain a precise measure of the transit depth. As a consequence, LD coefficients 
cannot be fitted, but are fixed to tabulated values obtained from stellar models 
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at different wavelengths and in the appropriate filters (e.g. Wakeford et al. 
and references therein). The James Webb Space Telescope will instead be placed 
at the Sun-Earth L2 point, and will not suffer from pointing interruptions. For 
this observatory, the wavelength dependence of the LD coefficients will be fully 
appreciated. 

Stellar activity features have been found to influence the measure of other 
transit parameters, too. A notable case is the one stellar density, which can be 
constrained to a high level of precision from planetary transits (Pont et al. 
Sozzetti et al. (2007). Studying the light curve of CoRoT-7b, where the individual 
transits are not fully resolved, Léger et al. discussed an underestimation 
of the stellar density from the transit fit compared to the derived spectroscopic 
value, likely due to stellar variability. Other studies, such as Alonso et al. (2009), 
Barros et al. One Oshagh et al. (2013), found effects on the apparent transit 
timing. Figure shows the simulated difference between a transit with and 
without starspot occultation, as a function of mid-spot time compared to mid- 
transit time, for the hot Jupiter WASP-10b (Barros et al. |2013). The largest 
effect is produced when the starspot bump is close to the transit edges, marked 
with vertical lines. In particular, variations in the mid-transit time (top panel) 
across multiple transits may mimic transit timing variations (TTVs), which might 
be misinterpreted as due to non-transiting companions (Agol et al. (2005). Barros 
et al. found hints of starspot contamination by observing a large number 
of transits with significant variations in the depth of individual transits. 

This latter is another indication that the study of planets around active host 
stars deserves a special attention. It requires a careful monitoring of the star to 
evaluate and correct for the contribution of the starspots to the stellar flux due to 
stellar rotation and their changing characteristics and location. 


7 Planetary transits and stellar granulation 


Stellar brightness inhomogeneities caused by stellar convection are another source 
of noise for transit parameter determination. Granulation produces variations on 
time scales from tens of minutes to hours. Hot plasma is uplifted towards the 
stellar surface and produces brighter areas, called granules by Dawes (1864) for 
the first time. Cooler plasma sinks between granules, and produces separating 
lanes between granules (Nordlund and Dravins Nordlund, Stein, and As- 
plund |2009). The first observations on granulation were carried out on the Sun 
eee but this phenomenon is now thought to be common among Sun- 


like stars (e.g. Chandrasekhar |1960} Schwarzschild Nordlund |1982; Dravins 


and Nordlund |1990; Gilliland et al. and refereces therein). 


7.1 Effect on measured stellar parameters 


The first direct effect of stellar granulation is on the stellar parameters, which 
are fundamental in order to extract the planetary radius from transit observables. 
The power of the granulation signal is correlated with the stellar surface gravity 
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Fig. 14. Simulated difference between a transit with and without starspot occultation, 
as a function of the mid-starspot time compared to mid-transit time, for the parameters 
of WASP-10b. From top to bottom: transit mid-time, orbital semi-major axis divided 
by stellar radius, planetary-to-stellar radius ratio and orbit inclination. From Barros 
et al. (2013), MNRAS 430, 3032-3047 (their Figure 10). 


g: as a star evolves from dwarf toward the giant phase, its g decreases; at the 
same time, its convection zones deepen and the granulation timescale increases. 
By measuring the photometric variations on a time scale of 8 hours, Bastien et 
al. highlighted brightness variations driven by granulation. The stars of 
this study span an effective temperature between 4500 to 6750 K, logg between 
2.5 and 4.5 (cgs units), and relative brightness variations lower than 1073. Taking 
advantage of g measurements obtained from Kepler asteroseismology observations 
(Chaplin et al. (2011), they derived the “flicker sequence” of stellar evolution. This 
sequence shows a tight correlation between the granulation-induced brightness 
variation and the stellar log g (see Figure 1 of Bastien et al. (2013). Within this 
phase space, Bastien et al. fitted an empirical relation (extended by Bastien 
et al.|2016) to infer g to within 25% from the granulation flicker of inactive Sun-like 
stars, from their main-sequence to their giant phase of evolution. 


In addition to surface gravity for the determination of stellar radii, stellar 
density is also particularly useful to correctly derive transit parameters. This could 
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be carried out via the so-called “asterodensity profiling” (Kipping [2014] Sliski and 
Kipping|2014). By using 588 catalogued Kepler target stars with asteroseismology 
measurements (Huber et al. Chaplin et al. (2014), Kipping et al. 
reported a linear trend between the stellar density and the 8 hr stellar flicker for 
stars with 4500K < Tog < 6500 K, 3.25 < logg < 4.43, and Kp magnitude < 14. 
The fitted relation has a model error in the stellar density of about 32%: this 
is an ~ 8x lower precision than asteroseismology, but the relationship can be 
applied to a ~ 40x wider sample of targets with no asteroseismic measurements. 
Capitalising on this, Bastien, Stassun, and Pepper re-estimated the stellar 
radii of 289 bright (Kp < 13) candidate planet host stars with 4500K < Tig < 
6650 K and assessed the impact on the parameters of their planets. They found the 
log g derived from broadband photometry and spectroscopy to be systematically 
lower that the one obtained from the flicker-based calibration (~ 0.2 dex median 
difference, with RMS = 0.3 for spectroscopy and 0.4 for broadband photometry). 
Hence, nearly 50 of the brightest stars in their sample resulted to be subgiants and 
their radii, together with those of their planets, on average 20%-30% larger than 
previously measured. 


7.2 Effect on transit parameters 


Stellar granulation affects the measured transit depth, duration determination of 
transit ingress and egress, and LD parameters. To quantify the effects, Chiavassa 
et al. (2017) used three-dimensional radiative-hydrodynamical stellar models from 
the STAGGER grid (Magic et al. to model the transit of a hot Jupiter, a hot 
Neptune and a terrestrial planet in front of the Sun and of a K star. 

The simulations were run by considering several convective turnovers and by 
taking into account temporal variations in the granulation intensity (e.g. about 
10 minutes for the Sun, Nesis et al. (2002). They were also performed in several 
wavelength windows, covering the bandpasses of most instruments used to ob- 
serve spectrophotometric transits. The modelled root mean square of the stellar 
flux in the visible (between 1 and 16 ppm) was close to the observed photometric 
variability of the SOHO quiet Sun data (10 to 50 ppm: Jenkins Frohlich 
et al.|2009) and was the largest in the visible. Then, different realisations of the 
stellar irradiance calculated from the granulation patterns were averaged in order 
to derive the photometric noise they produced. This resulted in detectable varia- 
tions of r,/R,, larger for the Sun than for the K star. In the visible, for example, 
they observed up to 0.90% and up to 0.45% uncertainties in the planet radius for 
the Sun and the K star simulations, respectively. Also, larger uncertainties were 
found for terrestrial than for giant planets. Figure illustrates the case of a 
terrestrial planet transiting in front of the Sun, in visible wavelengths, for differ- 
ent orbital inclinations. Given the required precision on planetary radii of a few 
percent at most, in order to be able to significantly constrain planetary core sizes 
(Wagner et al. (2011), such estimates imply that granulation must be considered 
a non-negligible source of uncertainty. Sulis et al. quantified the effect of 
granulation on the transit parameters, which can induce errors of up to 10% on the 
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Fig. 15. Simulated transit profile variations of a terrestrial planet transiting the Sun, due 
to different realisations of the modelled granulation pattern. Filled contours indicate the 
deepest and the dimmest transit profile for different orbital inclinations 7 in the 750-770 
nm band. Blue corresponds to i = 90.85°, green to i = 90.65°, yellow to i = 90.45°, and 
red to i = 90.25°. Credit: Chiavassa et al. (2017) (their Figure 11), A&A 597, A94, p. 
10, reproduced with permission @ESO. 


ratio between the planetary and stellar radius (r,/R,) for an Earth-sized planet 
orbiting a Sun-like star. Stellar flicker is already included in the error budget for 
the CHEOPS mission (Broeg et al. (2013), and it will be even more so for the 
PLATO mission (Rauer et al.|2014), which is designed to achieve a 2% precision 
on planetary radii. 


The apparent radius variation was also shown to be dependent on the orbital 
inclination of the planet, with additional effects on the duration of transit ingress 
and egress and correlations with the LD parameters. Morris et al. explored 
the correlations among these effects and found uncertainties on the planetary ra- 
dius up to 3.6% in the PLATO band, due to granulation and stellar oscillations 
only. These authors highlighted the importance of independent constraints on 
the transit impact parameter, or of follow-up observations at longer wavelengths, 
where LD is weaker, to reduce the significance of this effect. On the other hand, 
the method used to retrieve the transit parameters seems crucial in the determi- 
nation of the error budget: Gaussian processes, which are particularly suitable for 
the modelling of a stochastic phenomenon such as stellar granulation, showed to 
provide encouraging results with specific kernel choices (Barros et al. (2020). 
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In the context of the forthcoming James Webb Space Telescope (e.g. Gard- 
ner et al. Beichman et al. and the future Ariel mission (Tinetti et 
al. |2018), the potential impact of stellar granulation and oscillations on transit 
spectrophotometry has started to be evaluated. Sarkar et al. reported that 
the expected effect is weak, especially in the near-IR. There, its impact could even 
be negligible, but with a dependency on the atmospheric scale height and the 
transit duration. As a case which is most likely to be impacted by granulation, 
these authors mention the one of ‘a terrestrial planet with a secondary atmosphere 
orbiting a nearby Sun-like star on a long period’ (Sarkar et al. 2018, p. 2877). 


8 Conclusions 


In this review, we presented the main issues caused by stellar activity on exoplanet 
transit observations, as well as the main approaches which are used to correct for 
them. In several cases, the combination of ultra-precise measurements and theo- 
retical modelling is required to precisely quantify the possible biases on planetary 
parameters. Depending on the observations and the scientific goal, complementary 
observations can provide a useful way to better assess the effects of stellar activ- 
ity. On the other hand, attention to the determination of the physical properties 
of stellar activity features such as starspots, faculae and granulation and of their 
temporal evolution has grown in the past few years. The goal of a few percent 
precision in planet radii, and few tens of parts per million in transmission spec- 
troscopy, is challenged by the imprints of stellar activity on the observables, and 
requires these phenomena to be taken into account as an unavoidable component 
of star-planet observations. 

As we push our search towards terrestrial planets, the understanding of stellar 
activity and its proper correction becomes an increasingly necessary complement to 
improved instrument performance. The best outcomes will then result from syner- 
gies between different observation techniques and modelling perspectives. Thanks 
to these, we will be able to obtain as much precise as possible insights into the 
physical properties of the most interesting discovered systems. At the same time, 
this race towards small size planets with precise and accurate parameters will help 
refine our knowledge and understanding of stars. 


GB acknowledges support from CHEOPS ASI-INAF agreement n. 2019-29-HH.0. MD acknowl- 
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